Special scattering regimes for conical all-dielectric nanoparticles

All-dielectric nanophotonics opens a venue for a variety of novel phenomena and scattering regimes driven by unique optical effects in semiconductor and dielectric nanoresonators. Their peculiar optical signatures enabled by simultaneous electric and magnetic responses in the visible range pave a way for a plenty of new applications in nano-optics, biology, sensing, etc. In this work, we investigate fabrication-friendly truncated cone resonators and achieve several important scattering regimes due to the inherent property of cones—broken symmetry along the main axis without involving complex geometries or structured beams. We show this symmetry breaking to deliver various kinds of Kerker effects (generalized and transverse Kerker effects), non-scattering hybrid anapole regime (simultaneous anapole conditions for all the multipoles in a particle leading to the nearly full scattering suppression) and, vice versa, superscattering regime. Being governed by the same straightforward geometrical paradigm, discussed effects could greatly simplify the manufacturing process of photonic devices with different functionalities. Moreover, the additional degrees of freedom driven by the conicity open new horizons to tailor light-matter interactions at the nanoscale.

. Artistic representation of the considered silicon nanocone particles illuminated with a linearly polarized plane wave, showing the main advantage of the truncated cone geometry over the simpler shapes: the possibility of obtaining various important multipole effects within the same geometric shape.

Modes and multipoles
Multipolar expansion. Multipole decomposition is an important tool for analyzing the interaction of light with matter. Recently new expressions for multipole moments were introduced, valid for arbitrarily sized particles of any shape. Beyond the particular case of multipolar moments induced by an incident field in a structure, these expressions can be directly applied in the many areas, where the multipole decomposition of electric current density distributions are used 20 . Exact multipole moments (1)(2)(3)(4): where p α -electric dipole moment (ED), m α -magnetic dipole moment (MD), Q e αβ -electric quadrupole moment (EQ) and Q m αβ -magnetic quadrupole moment (MQ), α, β = x, y, z, J ω -current density inside the particle, j 1,2,3 -spherical Bessel functions of the first kind.
Hereinafter we will not consider any higher order terms 20 , since they are negligible in all the cases of interest.
In the far field region, the exact induced multipole moments produce an electric field given by 49 : Importantly, upon making the change n → −n in equation (5), the electric field produced by the ED does not change sign, while the electric field produced by the MD is reversed. For quadrupoles, the opposite happens, i.e., the MQ field is even while the EQ one is odd. The 'even-odd' character under space inversion is often referred to as the parity of the multipole field. This well-defined behavior under space inversion is at the core of all multipolar interference effects.
We also need to write down the following formula for the scattering cross-section (SCS) 20 , to directly compare the contributions of different multipoles to the SCS: Interestingly, we note that the SCS of each multipole is decoupled from the rest. Thus, every multipole can be viewed as an independent 'scattering channel' , through which the nanoparticle can exchange power with the environment 50 . Quasinormal modes. In the presence of an exciting field, dielectric nanoparticles have the ability to confine light within their volume. In other words, they act as subwavelength nanoresonators, supporting 'resonant modes' . This confinement becomes optimal near their resonant frequencies. To understand the nature of this phenomenon, one can trace an analogy with Fabry-Perot cavities. Once light enters the cavity, it bounces back and forth due to reflection from the walls until it finds a way to exit. However, for some special frequencies, a standing wave can form due to constructive interference, leading to resonant behavior. Another prominent analogy can be found in dielectric microcavities supporting whispering gallery modes, arising by total internal reflection after a roundtrip around the cavity. Resonances in dielectric nanoparticles can be visualized as a mixture between the two examples above; they can arise from standing waves along their transverse dimension (similar to whispering gallery modes), their longitudinal one (resembling Fabry Perot cavities), or a combination of both.
In an ideal system, such as a Fabry Perot resonator bounded by two perfectly conducting mirrors, light can never escape and will oscillate back and forth from the walls for an infinite time. The resonant frequencies in this context are real numbers. However, in an actual physical system, energy eventually dissipates into the environment, either by absorption or radiation. The resonant frequencies are then complex, implying that no stationary excitation can fully access them. Despite this, the system's response (e.g., scattering) is strongly affected when approaching a resonant frequency along the real axis. Thus, understanding their behavior and radiation characteristics is of great importance for the design of nanophotonic devices. Formally, resonant modes are eigensolutions of Maxwell's equations, supplemented with a set of boundary conditions. For an isolated dielectric (1) p α = − 1 iω d 3 rJ ω α j 0 (kr) + k 2 2 d 3 r 3(r · J ω )r α − r 2 J ω α j 2 (kr) (kr) 2 , Q e αβ = − 3 iω { d 3 r 3(r β J α + r α J β ) − 2(r · J ω )δ αβ j 1 (kr) kr + +2k 2 d 3 r 5r α r β (r · J ω ) − (r α J β + r β J α )r 2 − r 2 (r · J ω )δ αβ j 3 (kr) (kr) 3 }, (4) Q m αβ = 15 d 3 r r α (r × J ω ) β + r β (r × J ω ) α j 2 (kr) (kr) 2 , www.nature.com/scientificreports/ nanoparticle in a homogeneous environment, the boundary conditions can be replaced by a radiation condition. Assuming a dispersionless permittivity ε(r) , we write Maxwell's equations with a time dependence of the form e −iωt (we omit for brevity the spatial dependences in the fields and the permittivity): together with the radiation condition: where r is a unit radial vector and ω m = ω m − iγ m is the resonant frequency or eigenfrequency associated with the eigenmode E m H m . Equation (8) essentially tells that the mode field must behave as a spherical wave at infinity, i.e. follow a dependence E m ∼ E 0 r e iω m r/c /r . Interestingly, to ensure energy decay in the time domain, γ m is always a real, positive number. Therefore, far from the origin, the mode amplitude blows up due to a term e γ m r/c /r . In consequence, modes in open systems cannot be normalized by standard means, and often receive the name 'quasinormal' modes (QNMs). Despite this, the scattered field everywhere inside the nanoparticle can be accurately described as a linear combination of the contributions of several QNMs, i.e.
where α m (ω) describes the coupling of the m-th QNM to the incident field 51 . Multipolar resonances in dielectric nanoparticles are associated with a QNM. In the case of spheres, every resonant mode radiates as a specific multipole. However, this is no longer true for arbitrary shapes. In a general setting, the scattering pattern of a QNM can be described as a mixture of multipole contributions. To determine how much a QNM 'matches' a given multipole, (i.e. its multipolar content) one can define the m-th 'eigen' current J m = −iω m ε 0 ε p − 1 E m ,where ε p is the relative permittivity of the nanoparticle. Introducing J m into equations (1)(2)(3)(4), and evaluating them at the complex frequency ω m yields the desired multipolar content of the m-th QNM 52 . At a real frequency ω , the QNM will radiate as a combination of its intrinsic multipole moments.
However, the particle symmetry imposes strict bounds on the multipolar contents of its QNMs. For instance, cylindrical symmetry prevents a mixture of multipoles with even and odd parity (refer to Fig. 2a-c). This has important implications, e.g., for the design of Kerker meta-atoms: since the electric and magnetic dipoles have opposite symmetry, the radiation pattern by a single mode cannot be directional. Thus, engineering a resonant Kerker effect in a dielectric nanodisk requires overlapping the resonances of two QNMs with electric and magnetic dipolar character, respectively.
Unlike cylinders, truncated nanocones lack the vertical mirror plane, as schematically depicted in Fig. 2d. As a result, the QNMs can radiate as a combination of multipoles with even and odd parity, such as the electric www.nature.com/scientificreports/ and magnetic dipole. Consequently, the directivity of a single QNM can be enhanced in the forward or the backward direction, as can be appreciated from the radiation patterns of the two lowest order QNMs of a nanocone ( Fig. 2e-f). This also immediately implies the appearance of a bianisotropic response. Thus, truncated nanocones not only hold an additional degree of freedom (conicity) with respect to more studied geometries such as cylinders: conicity provides a simple strategy to control the even-odd mixture of multipoles in a QNM, leading to new exotic effects, such as single mode directivity, strong coupling and exceptional points 53 .

Kerker effects
In this section, we talk about various types of Kerker effects (generalized and transverse Kerker effects) in silicon truncated nanocones. The Kerker effect became more commonly used recently due to the explosive growth of dielectric nanophotonics. In this regard, it becomes necessary to study this type of effects considering different scatterer geometries. The Kerker effect is a unidirectional forward or backward scattering (Generalized Kerker effects) or enhanced side scattering (transverse Kerker effect). The multipole decomposition can serve as an excellent tool to study the Kerker effect. It becomes possible to show the electric field as the sum of multipole contributions, where each term is responsible for the field of a particular multipole 20 . Generalized Kerker conditions. For the first time, the Kerker effect was discovered for spherical particles with the dielectric constant equal to the magnetic permeability. In such particles, when the electric field amplitudes of the electric and magnetic dipoles are equal, as well as at a certain phase difference of the dipoles, one can observe only forward or backward scattering 21 . There are many works that show the experimental realization of such effect in dielectric nanoparticles [54][55][56] . Over time, other meaningful combinations of multipoles and the phase difference between them were discovered, so the Kerker effect had to be "generalized" 25,57 . Nowadays, the term "generalized Kerker effect" is used, when it comes to pronounced forward or backward scattering.
Amplitude ratio for generalized Kerker conditions. To describe forward or backward scattering we can write the scattered field of an arbitrary shaped particle under x-polarized light 31 , inserting n = (0, 0, n z ) in (5), we obtain (9): The direction of forward or backward scattering can be linked to the unit vector n. Let us set n = (0, 0, 1) for the forward scattering, and n = (0, 0, − 1) for the backward. We consider the case of a plane wave illumination with x-polarization.
Let us assume that we do not have backscattering, and all the multipole contributions that do not participate in the Kerker effect tend to zero. In such cases, the following expressions are valid: It should be emphasized that the last combination (14) includes all multipoles up to quadrupoles, and therefore is a superposition of the previous cases. Also, it is worth noting that combination of four multipoles can lead not only to the generalized Kerker effect, but also to the transverse Kerker effect. Such multipole combinations will be discussed in Sect. "Conclusions on the Kerker effects".
Knowing the multipoles values for 'no backscattering' cases, we can express the forward scattering amplitudes: www.nature.com/scientificreports/ Further, considering the equations (12,(4)(5)(6)(7)(8)(9)(10)(11), we obtain the ratios between the multipole contributions for each combination of moments leading to the Kerker effect (see Supplementary 2): Phase difference for generalized Kerker conditions. Another important Kerker effect condition is the phase difference between the interacting multipoles. For the original problem of Mie scattering on spherical particles, it was shown that electric and magnetic multipoles of the same order have opposite parity with respect to cos θ (θ is the scattering angle with respect to the forward direction) 58 . This means that simultaneous forward scattering amplification and backscattering suppression can be achieved not only for ED and MD interference, but also for any higher-order multipoles, and scattering can be further enhanced or suppressed 59 .
The fact is that the phase will directly affect the shape of the scattering diagram we obtain. For all the above combinations of multipoles, specified phase differences can be found to ensure forward scattering.
In our work, all the generalized Kerker conditions with forward direction were obtained for truncated conical nanoscatterers (Fig. 3). Figure 3 shows the Kerker effects calculated numerically, using Comsol Multiphysics, and analytically, using equation (5). With the help of this equation, it is possible to construct the far field distribution by varying the values of the multipole moments in such a way as to obtain the Kerker effect. The dependence of the multipole contributions to the scattering cross-section on the radiation wavelength is shown in Fig. 3(a-e). The points of the Kerker effect are marked in red, following the equations (19)(20)(21)(22). To confirm this, the analytical (f-j) and numerical (k-o) far-field patterns are shown for each case. The scattering shape is similar but not completely the same due to the minor contributions of other multipoles in the numerical calculations.
Thus, various kinds of Kerker effects, and necessary far field and scattering cross-section conditions are shown for real silicon truncated cones. The resulting shapes of radiation patterns for cones, obtained in COMSOL Multiphysics, are in a good agreement with the expected Kerker-type shapes for an ideal point calculations. In real systems, it is impossible to obtain the ideal combination of the necessary multipoles, and therefore weak discrepancies are observed. The absence of the need to use higher-order multipoles is also evidenced by good coincidence of the sum of scattering cross-section of multipoles with calculation of scattering cross-section, using the integration of the Pointing vector over a closed surface in the far-field zone and the normalization to the incident field intensity ( Fig. 3a-e, gray and orange lines).

Transverse Kerker conditions.
Another equally important scattering feature is the transverse Kerker effect, which was first described in 29,60 . The main cause of the transverse Kerker effect is the combination of multipoles featuring scattering to the sides only leaving the small portion of forward scattering according to optical theorem 61 . Such effect can be obtained both through simple combinations of two multipoles 62 and through more complex configurations 30 . www.nature.com/scientificreports/ Let us assume that both backward and forward scattering are suppressed, and all the multipole contributions not participating in the transverse Kerker effect tend to zero. Then, the following system of equations is valid for possible ways to obtain transverse Kerker effect (4 multipoles and 2 multipoles): After some simple algebra, one can obtain the following conditions for the multipole amplitudes: For all the above combinations, phase differences between multipoles can be obtained. Transverse phase conditions are: In our work, all the transverse Kerker conditions are shown for truncated conical nano-scatterers (Fig. 4). Figure 4 shows the different transverse Kerker effects obtained both numerically and analytically, calculated in the same way as in Fig. 3. Scattering patterns shape is an obvious reason why this effect is called "transverse". By the proper choice of the multipole phases and amplitudes, it is possible to achieve side scattering together with forward and backward scattering suppression.
Conclusions on the Kerker effects. Kerker effects are indeed unique phenomena, which are key to a range of novel photonic devices. Prior to our work, each one of these had previously been only shown for completely different shapes and materials of nano-scatterers. The following are some examples demonstrating various Kerker-type effects (Table 1).
In this section, we have demonstrated for the first time all known Kerker effects for single nanoscatterers within the same geometry. Thus, nanocones represent a versatile, fabrication-friendly platform for the implementation of new photonic devices benefitting from a comprehensive toolbox of multipolar interference effects including flexible tailoring of scattering patterns.

Non-scattering regimes: anapole and hybrid anapole
In the past few years, the emergent field of 'anapole electrodynamics' is experiencing exponential growth 69 . Anapoles are semi-nonradiating sources that arise due to the destructive interference of the quasistatic electric dipole moment and the toroidal dipole in the far field. Alternatively, in a more general picture, they can be understood as being originated by the destructive interference of symmetry-compatible quasinormal modes 33 . However, the energy stored by the quasinormal modes within the nanoparticle is nonzero, leading to counterintuitive light-matter interaction processes in the absence of elastic scattering. Until very recently, the suppression of scattering in these states was limited to the electric dipole contribution to radiation. The experimental demonstration of hybrid anapoles (HA) following the pioneering theoretical proposal 70 , has evidenced the possibility to simultaneously overlap the zeros of all the dominant multipolar channels through a careful design of the nanoparticle geometry. These novel states are much more promising than their dipolar counterparts for a number of reasons 33,34 ; despite the larger volume of the nanoresonator required to obtain them, the scattering suppression is improved by more than 20 times, while the excited quasinormal modes store approximately 10 times more energy. Such values exceed by far the performance of anapoles and 2 nd order anapoles 71 arising in homogeneous disk nanoresonators.
In nanocylinders pertaining to the cylindrical symmetry, it was shown in 33 how the anapoles from different multipolar orders but equal parity were connected. Counterintuitively, such connection allows to simultaneously www.nature.com/scientificreports/ overlap four anapoles with just the two geometrical degrees of freedom of the nanoparticle. As mentioned earlier, however, the even and odd multipoles are no longer necessarily coupled once the reflection symmetry in z is broken. In this section, we aim at investigating the effect of z-symmetry breaking on a HA nanocylinder, by introducing a small geometrical perturbation on the upper radius. Such a situation occurs quite often, (unintentionally) during a sample fabrication process. It is rarely possible to obtain particles of an ideal shape, and samples with small defects can often be found. In practice, the manufactured nanocylinders are most likely truncated cones 72,73 . Figure 5a displays the exact multipole decomposition of the HA regimes of the silicon cylinder studied in 33,34 . The simultaneous suppression of the ED, MD, EQ, and MQ channels at the same spectral position can be appreciated. The inset depicts the norm of the electric field at the spectral point with the lowest scattering, hereafter referred to as the 'HA wavelength (λ HA )' . The field can be seen to be strongly concentrated within the nanoparticle, (where two hotspots appear), with the exception of a few hotspots at the surface with a size of the order of λ HA /40. Thus, even in the near field outside the particle, the incident wave is barely distorted by scattering. We now keep the height of the resonator constant and progressively increase conicity ( 1 − R top /R bottom ). The multipole decompositions for three selected conicities are shown in Fig. 5a-c In all cases, illumination from the top is considered. A log scale is used to enhance the contrast of the zeros in the spectrum.
Firstly, we observe a spectral blue-shift of all resonant features. This is expected, due to the overall size reduction of the lateral dimensions of the resonator. Secondly, the original MQ and ED anapoles blue-shift faster  www.nature.com/scientificreports/ than the EQ and MD. As a result, the four anapoles are no longer superposed in the spectrum, and the scatterer becomes less 'transparent' to the incident illumination. This is demonstrated by calculations of the total SCS as a function of conicity (Fig. 5d). However, we also observe a progressive emergence of new anapoles in the ED and the MQ SCSs (red and purple solid lines, respectively), coinciding with the maxima in the EQ and MD SCSs near the HA. These new features become more pronounced with increasing conicity. Interestingly, their formation is acompanied by a 'split' (two separated anapoles close to each other, see Fig. 5c) in the unperturbed ED and MQ resonances.
We evaluate quantitatively this effect by calculating the wavelength shifts of the multipolar anapoles as a function of conicity (Fig. 6). The white-dashed lines show the paths followed in parameter space (only the anapoles involved in the vicinity of the HA regime are investigated). The unperturbed anapoles in the ED and MQ channels are labeled as A 1 and B 1 , respectively (Fig. 6a, d). As could be anticipated from the results in Ref. 33 , the latter follow the same path in parameter space. This is because they are both associated to the same resonant QNM. The same occurs with the MD and EQ anapoles (Fig. 6b, c).
For values of 1 − R top /R bottom > 0 , we confirm the appearance of new anapoles in both the ED and MQ channels, manifesting as pronounced dips in their contributions to the SCS (Fig. 6a, d). For convenience, we denote them as A 2 and B 2 . Importantly, we notice that they both follow a similar path in parameter space. What's more, there seems to be a connection between the paths of A 2 and B 2 and those followed by the original EQ and MD anapoles; after their appearance, the former spectrally overlap with the latter. In contrast, no appreciable new features can be seen when 1 − R top /R bottom < 0 . In a first regard, this observation seems contradictory, since symmetry is analogically broken. However, later the physical reason for it will become clear. It is necessary to mention that in the calculations of the conicity we keep R bottom constant and increase R top , then the volume of the 1 − R top /R bottom < 0 is larger than 1 − R top /R bottom > 0.
To gain more insight, we briefly revisit the origin of anapoles, from the QNM perspective 74,75 . As an example, consider a scatterer supporting one resonant QNM radiating to the ED channel. The SCS would then be well approximated by the first term in equation (6) www.nature.com/scientificreports/ Here, p 1 is the induced ED moment by the resonant QNM, and p bg corresponds to a non-resonant contribution stemming from QNMs outside the spectral range of interest, as well as direct scattering from the object's shape 76 . A zero (anapole), takes place when p 2 = 0 , so that the resonant contribution cancels out with the background: We note that, in general, p bg is a smooth function of frequency that cannot, a priori, be controlled by design, while p 1 becomes non-negligible only near the resonance frequency. In particular, one can model p 1 as a Lorentzian function centered at the resonance frequency. In the presence of a second resonant QNM, equation (40) can be expressed as: Equation (41) provides one more degree of freedom to reach the anapole condition. In principle, it is possible to tune p 1 so that p 1 = −p bg − p 2 , or change p 2 so that p 2 = −p bg − p 1 . Since p 1,2 vary naturally in amplitude and phase near the resonance frequencies of the associated QNMs, a new anapole should emerge for every new QNM contributing to the ED cross-section 77 . This conclusion can be readily generalized to any multipole.
From the above discussion, it follows that anapole formation is intrinsically related to modal evolution. In particular, the emergence of a new anapole in the ED and the MQ SCSs must be directly connected to an additional QNM contribution to these multipoles. With this idea in mind, we calculate the resonant QNMs in the vicinity of the HA as a function of the conicity (Fig. 7a). Three QNMs were considered: the two QNMs primarily responsible for the HA effect (M1 and M2 in Fig. 7a) and, for completeness, a third QNM lying within the visible range (M3 in Fig. 7a).
In Ref. 33 , it was shown that M1 was associated with the formation of both the ED and MQ anapoles ( A 1 and B 1 ), while M2 was responsible for the MD and EQ ones. This is further confirmed here by comparing the paths of the anapoles in Fig. 6 with the evolution of the resonant wavelengths in Fig. 7a. The role of M3 can be safely neglected in the discussion (or considered part of the background), since it is spectrally isolated from M1 and M2 in the range of parameters considered.  www.nature.com/scientificreports/ We now focus our attention on the multipolar character of QNM M2 (Fig. 7b). We calculate the intrinsic multipole moments as explained in Sect. "Quasinormal modes". This indicates us whether a QNM will contribute or not to the SCS of a given multipole. For our purposes, if the contribution of a specific QNM is zero, it cannot play a role in the formation of an anapole in that multipole SCS.
The results reveal a surprising fact: for an unperturbed cylinder, the contributions of M2 to the ED and MQ are zero, while the MD and EQ are not, as expected since it interferes in those two channels to form anapoles. However, when 1 − R top /R bottom � = 0 , M2 starts contributing to the ED and MQ SCSs. The MQ and ED contents of M2 grow as a function of the perturbation. It immediately follows from the discussion above that a new anapole can emerge in the two multipoles under consideration. Furthermore, the contributions to the MD and EQ SCS also change: for 1 − R top /R bottom < 0 the QNM is better matched to the EQ. Conversely, for 1 − R top /R bottom > 0 , the MD becomes dominant.
The fact that M2 couples to multipoles with opposite parities when symmetry is broken can be understood from Fig. 2 and the discussion in Sect. "Quasinormal modes". In a truncated cone, the modes are no longer eigenstates of parity, and therefore can scatter light as a mixture of even and odd multipoles.
In brief, the emergence of anapoles A 2 and B 2 is due to mode M2 being able to contribute to the ED and MQ SCS, once cylindrical symmetry is broken. This also explains the apparent connection between the paths followed by the EQ and MD anapoles (associated with M2) and the new emerging anapoles. The reason why A 2 and B 2 appear only for 1 − R top /R bottom > 0 can be elucidated from the behavior of the resonant wavelengths in Fig. 7a. It can be seen that M1 and M2 remain spectrally close when increasing 1 − R top /R bottom . In contrast, for the range of parameters considered, M1 rapidly redshifts towards the near IR with decreasing 1 − R top /R bottom , while M2 remains in the visible. In consequence, M1 is spectrally isolated from M2, and the latter does not play an important role anymore. Thus, no new anapoles are formed.
Regarding the robustness of the HA regime under small perturbations, extensive numerical tests demonstrate that it is still possible to recover the spectral overlap of the four anapoles up to perturbations in the range of 1 − R top /R bottom = {0.2 , −0.2} by varying the resonator height (not shown). Interestingly, the lower limit appears to coincide with the region where A 2 and B 2 become more pronounced (Fig. 7a,b). This suggests that, once the MQ and ED content of M2 becomes significant, the MQ and ED anapoles can no longer be tuned independently from the EQ and MD ones, which constituted the general strategy to design HA in nanocylinders.
Next, we study the possibility to realize HA in truncated cones for conicity close to 1 (Fig. 8). Keeping the bottom radius and the height constant, we perform calculations of the total SCS for conicities in the range between 1 (perfect cone) and 0.5 (Fig. 8a). Our simulations reveal the existence of regions with strongly suppressed SCS (white square in Fig. 8a). Figure 8b displays the multipole decomposition for a nanocone with conicity 0.88. A HA can be seen to form due to the overlap between the ED and EQ anapoles, and the spectral proximity of MD and MQ anapoles. For the example at hand, scattering is suppressed by approximately 8 times the average SCS in the visible range. In the wavelength featuring the lowest SCS the scatterer is virtually transparent to the incident plane wave, as demonstrated in the field profile shown in Fig. 8c. The incident plane wave is seen to propagate undistorted by the cone. Although significant, the result leaves room for improvement. More intensive multiparameter searches can lead to a closer overlap between the four anapoles, and an even more pronounced scattering suppression.
The study in this section reveals the strong effect that small perturbations in conicity can have in the resonances responsible for the HA, since they effectively break the underlying symmetry, and render previously closed channels open, which leads to additional interactions between the modes. The knowledge of these can be important for future applications that benefit from the exotic properties of HA, for instance, metasurface engineering, sensing, etc. Furthermore, we have demonstrated, for the first time, the possibility to achieve HA in www.nature.com/scientificreports/ nanocones. In this way, we have expanded the library of all-dielectric nanostructures that support HA, beyond cylinders 33 and ellipsoids 36 . Once again, cones present themselves as a simple, fabrication-friendly platform for the implementation of multipolar interference effects.

Superscattering regime
In the last section, we introduce the superscattering effect. In contrast to the hybrid anapole this effect concerns scattering enhancement from a subwavelength particle. This unique regime has already found a plethora of emerging applications, such as, e.g. sensing 78 , energy harvesting 79 , radar deception 80 , etc. To enhance scattering, it is generally believed that one needs to spectrally overlap the resonant frequencies of several QNMs scattering to different multipoles 37 . It follows from a simple geometry tuning where the thicknesses of core-shell structures were optimized to precisely bring modes together. Since the structure retains spherical symmetry, no multipole mixing is allowed, and thus there is no interaction between the different QNMs. Several other works followed on the same foot-steps; the superscattering regime has been investigated in a variety of structures retaining spherical symmetry 39,58,[81][82][83] and was recently confirmed experimentally in the microwave frequency range 84 .
The generally accepted definition of superscattering is as follows 37 : the total scattering cross-section of a resonator must exceed (by far) the maximum scattering cross-section of a dipolar particle with spherical symmetry 37 . The maximum contribution of a multipolar channel to the scattering for the specified resonator can be formulated as 37,85 , where l is the total angular momentum and λ is the wavelength. Thus, for the first order multipoles (electric and magnetic dipoles) the maximum cross-section reduces to Hereafter we refer to this limitation as the dipolar maximum for the spherically symmetric scatterer or the DM as an abbreviation. Hence, superscattering corresponds to the condition C max >> C p max . Recently, it has been noticed that departing from the spherical scatterer to a scatterer lacking a rotation and/ or a reflection symmetry, this bound no longer holds 85 . Therefore, the scattering maxima of the multipoles are no longer limited by equation (42). In fact, to our knowledge, in a scatterer of arbitrary size and shape, the contribution of a multipole to the scattering cross-section is not bounded 50 . In the literature, however, the scattering limitations defined for spherical scatterers remain as a benchmark for the superscattering occurrence. Henceforth, we are normalizing the scattering by the dipolar channel maxima defined in equation (43) to establish consistent comparisons for the scattering enhancement in different setups.
In our recent work 85 , we have suggested and experimentally demonstrated a new paradigm to achieve superscattering. It benefits from the vastly growing fields of non-Hermitian physics. Indeed, it has been established that subwavelength open-cavity resonators support QNMs that may interact and collectively lead to a plethora of light-matter interactions effects such as enhanced directionality or broadband scattering.
We have shown 85 starting from a spherical resonator supporting two QNMs that radiate independently to two multipolar channels of different order but equal parity, a fine tuning of the vertical and/or the horizontal radius (i.e., deforming the sphere into a spheroid) results in increasing scattering on a certain multipole so that it exceeds the limitation introduced in equation (42). It also proves that breaking the fully symmetric structure, in this case defined by the sphere, into an axisymmetric structure such as the ellipsoid may lead to strong coupling www.nature.com/scientificreports/ between two QNMs and furthermore allows controlling the multipolar character of the scattered wave (i.e. its far field characteristics). Thus, for the first time the concept of super multipoles was introduced. This newly suggested pathway to superscattering opens uncharted territories for designing resonators that capture incident photons in an unusually large area. Compared to core-shell spherical particles where superscattering can be achieved only through the accidental overlapping of modes, scattering of reduced symmetry resonators can be enhanced both accidentally (spectrally tuning the resonance frequencies of two noninteracting QNMs) or by carefully tailoring the interference of coupled QNMs. As emphasized in the previous sections, the tuning of the conicity offers more degrees of freedom in controlling the radiation patterns of the QNMs that are unavailable for scatterers with cylindrical symmetry. Lack of the reflection symmetry in conical geometries allows for interference between modes with unequal parities. This additional degree of freedom has an important consequence on realizing super multipoles as we later explain. Figure 9a displays the total scattering cross-section normalized by the DM (Eq. (43)) as a function of the conicity and the wavelength. To reduce the computational complexity, we again assume dispersionless Si materials similarly to the setup described in the previous section. Cone's geometrical parameters are enlisted in the figure caption, we here selected a subwavelength cone in the visible and near-infrared region as a practical low-loss superscatterer structure. It is apparent that the normalized scattering exceeds by more than double the DM in the specified range. One can also observe that the now presented "superscattering" broadens while tuning conicity with multiple appreciable superscattering peaks. Noteworthily, the normalization factor (the DM) is not a function of the geometry, therefore, one can assume that the value of superscattering can be readily increased by considering factors of geometry along with the operating wavelength range and the material composition of the resonator.  Fig. 9b show the internal field distribution of the two modes, (detailed field profiles inside and in the vicinity of the cylinder are shown in supplementary figure S2). M 1 ' and M 2 ' avoid the crossing, and redshift as the cone volume increases (for an increased conicity). We consider these two QNMs as the focus of our investigation to alter and potentially enhance scattering.
From Fig. 9d, it can be appreciated that the two even multipolar channels, the electric dipole and the magnetic quadrupole, contribute to the scattering at spectral points coinciding with the resonant wavelength of M 2 ' . We notice that the electric dipole is showing superscattering where it exceeds the DM by about 1.5 times. The latter are closely related to the super dipole resonances studied in Ref 85 . Similarly, M 1 ' and M 3 ' appear to radiate on the odd multipolar channels, the magnetic dipole and the electric quadrupole, respectively. The magnetic dipole displays a scattering cross-section with values slightly above the DM. Summarizing, mode M 2 ' is even and scatters as a combination of electric dipole and magnetic quadrupole moments, while odd QNMs M 1 ' and M 3 ' are odd and scatter as a combination of magnetic dipole and electric quadrupole moments.
We show in Fig. 9c and d cases for conicity 0.25 and − 0.4, respectively. In the first, the electric dipole crosssection is now almost double the DM, while the magnetic dipole is reduced below that level. When decreasing the upper radius of the cavity, the resonant wavelengths of QNMs M 1 ' and M 2 ' are blueshifted [ Fig. 9b]. As a result, the multipole resonances are also blueshifted [ Fig. 9c]. The opposite occurs when increasing the upper radius, [Fig. 9e]. However, a second magnetic dipole scattering peak is now observable, overlapping with the electric dipole and magnetic quadrupole peaks. Meanwhile, we show for all the three cases the scattering patterns at the distinguished dipolar scattering peaks. The scattering is more directive when both dipoles are present and comparable in values, and it takes obviously a dipolar scattering signature when either of the dipoles is dominant. To have a clearer picture of the mechanism at play that controls the position and amplitudes of the superscattering multipolar channels, we show in Fig. 10 the scattering evolution of the first four multipoles as a function of the conicity and the wavelength.
When compared to the QNMs' dispersion in Fig. 9b, we notice that the power scattered by the electric dipole is mainly attributed to mode M 2 ' while it appears that modes M 1 '-M 3 ' are able to scatter through the magnetic dipole channel. For this reason, we see that the magnetic dipole appears with two peaks on the scattering spectrum where it is especially noticeable for negative higher conicities (conicity < − 0.2). Modes M 3 ' and M 4 ' across www.nature.com/scientificreports/ the parametric space are radiating on the electric quadrupole channel. One can notice also that the magnetic quadrupole channel is mostly apparent aligning with the electric dipole peak scattering. Although the quadrupolar channels are not superscattering within themselves in the current design, they have important contributions to the total scattering cross-section. It can be concluded that modes M 1 ' and M 2 ' of the cylinder are allowed to interfere with each other after breaking the vertical mirror symmetry i.e., with conicity ≠ 0. This implies that both modes are now able to radiate on the even and odd parity multipolar channels. One important result of this interference is an increase of the electric dipolar superscattering with almost doubling the DM at small conicities. We can also appreciate that M 2 ' start to scatter on the magnetic dipolar channel in the bigger conicities region of the parametric figure. Considering this mode is originally even parity in the cylinder, it becomes hybrid when tuning conicity in the sense that it scatters as a mixture of electric and magnetic dipole moments. Therefore, a single QNM in the conical geometry is able to drive opposite parity multipolar channels to superscattering regime. This type of superscattering cannot be obtained in cylindrical resonators without spectrally overlapping the resonant wavelengths of at least two QNMs. Such hybridization of superscattering channels is particularly important for the superscattering effect since it also decisively impacts the directivity of the overall scattering. This means, that it can be utilized to enhance forward or the backward scattering upon prerequisites. Another important consequence of M 1 '/M 2 ' interference is the broadening of the superscattering. This can be seen in both the lower and upper regions of the interference.
The above systematic exploration of the superscattering achievable through designing conical geometry has revealed a new approach to obtain desired outcomes in superscattering features. We focused on tuning two modes (M 1 ' and M 2 ') of a cylinder that are in spectral proximity, sufficiently confined within the resonator domain, and with an unequal parity. This strategy can be sufficient but is not necessary to achieve superscattering in the geometry under consideration. However, the investigated mechanism highlights important features of superscattering in a cone, namely, (i) the possibility to exceed the DM with the help of the novel supermultipoles and (ii) the mixed character of the involved QNMs, which scatter as a combination of even and odd multipoles, unlike in spheres or cylinders. Alternatively, other strategies remain possible as well, however with optimized conical structure in this work we believe that several important aspects of the superscattering conical geometry have been revealed.
Thus far we have shown a remarkable control over scattering, due to the design flexibility and symmetry rules governing the conical geometry. Tuning a single geometrical parameter (in this case-the conicity) can result in strong coupling between two QNMs of opposite parity, as was previously shown very recently 53 . Remarkably, this mechanism can be utilized to enhance scattering and manipulate the radiation pattern at the resonance to render it more directive. It indeed provides a plethora of opportunities to optimize truncated cones for applications in optics.

Conclusion
In this work, we have proposed the conical geometry as a universal platform for obtaining many important optical effects that are of interest in nanophotonics. For the first time, all known Kerker effects were obtained on one scatterer shape for a real dispersion of the refractive index of silicon. It was found that the recently discovered non-scattering hybrid anapole regime can also be obtained in conical nanoparticles, in addition to elliptical and cylindrical ones. It is shown how, by changing the geometric parameters, it is possible to adjust the anapole regime of different multipoles and obtain a hybrid anapole regime. We also studied the possibility of obtaining superscattering, the effect inverse to the anapole regime, in conical particles, showing the dependence on the geometry of the scatterer. This research takes a step towards nanophotonics of more complex shapes with the ability to fine-tune the effects that can be obtained on a single shape of a nano-scatterer. This research significantly reduces the cost of developing photonic devices and opens up new horizons for the practical application of the next-generation photonics. This work can find its application in various fields of research, for example, to create various dielectric nanoantennas in the form of a chain of resonators in comparison with chains of spheres or cylinders, or for metasurfaces based on which it will be possible to obtain a number of optical effects that were previously inaccessible.

Data availability
All data generated or analysed during this study are included in this published article [and its supplementary information files].